Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 November 201 1 (MN KTbX style file v2.2) 



Gamma Ray Bursts in the comoving frame 



G. Ghirlanda ^*, L. Nava^, G. Ghisellini\ A. Celotti^, D. Burlon^ S. Covino\ A. Melandri^ 



^INAF— Osservatorio Astronomico di Brera, Via E. Bianchi 46, 1-23807 Merate, Italy 

^SISSA - via Bonomea, 265, 1-34136 Trieste, Italy 

^ Max-Planck-Institut fUr Extraterrestrische Physik, Giessenbachstrafie 1, D-8547S Garching, Germany 



o 

(N 

-(— > 
o 

O 
m 

m 

X 
6 

> 

O 
O 



% 



ABSTRACT 

We estimate the bulk Lorentz factor Fq of 31 GRBs using the measured peak time of their 
afterglow light curves. We consider two possible scenarios for the estimate of To: the case of 
a homogeneous circumburst medium or a wind density profile. The values of Fq are broadly 
distributed between few tens and several hundreds with average values ^-^138 and ^66 for 
the homogeneous and wind density profile, respectively. We find that the isotropic energy and 
luminosity correlate in a similar way with Fq, i.e. i^isoC)cFo^ and LisoOcFo^, while the peak 
energy i^poakOcFo. These correlations are less scattered in the wind density profile than in the 
homogeneous case. We then study the energetics, luminosities and spectral properties of our 
bursts in their comoving frame. The distribution of L[g^ is very narrow with a dispersion of 
less than a decade in the wind case, clustering around L[^^^ 5 x 10''* erg s^^. Peak photon 
energies cluster around E' -^^ 6 keV. The newly found correlations involving Fq offer a 
general interpretation scheme for the spectral-energy correlation of GRBs. The i^pcak — E\so 
and i^pcak — Liso correlations are due to the different Fq factors and the coUimation-corrected 
correlation, iJpcak — E-^ (obtained by correcting the isotropic quantities for the jet opening 
angle 6'j), can be explained if 6'?Fo = constant. Assuming the E'pcak — E-^ correlation as valid, 
we find a typical value of OjVo^ 6-20, in agreement with the predictions of magnetically 
accelerated jet models. 

Key words: Gamma-ray: bursts — Radiation mechanisms: non thermal 



1 INTRODUCTION 



The discovery of the afterglows of Gamma Ray Bursts (GRBs - 
Costa et al. 1997) allowed to pinpoint their position in the X-ray 
and Optical bands. This opened a new era focused at measuring the 
spectroscopic redshifts of these sources. The presents collection of 
GRBs with measured z consists of 232 events. In 132 bursts of 
this sample (updated in this paper) the peak energy Sp^J^i^ of their 
vF^ prompt emission 7-ray spectrum could be constrained. In turn, 
for these bursts it was possible to calculate the isotropic equivalent 
energy Eiso and luminosity Liso- The knowledge of the redshifts 
showed that two strong correlations exist between the rest frame 
peak energy -Epcak and Eiso or Liso (also known as the "Amati" 
and "Yonetoku" correlations - Amati et al. 2002, Yonetoku et al. 
2004, respectively). 

The reality of these correlations has been widely discussed in 
the literature. Some authors pointed out that they can be the result 
of observational selection effects (Nakar & Piran 2005; Band & 
Preece 2005; Butler et al. 2007, Butler, Kocevski & Bloom 2009; 
Shahmoradi & Nemiroff 2011) but counter-arguments have been 
put forward arguing that selection effects, even if surely present. 
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play a marginal role (Ghirlanda et al. 2005, Bosnjak et al. 2008, 
Ghirlanda et al. 2008; Nava et al., 2008; Krimm et al. 2009; Am- 
ati et al. 2009). The finding that a correlation Ep(t)-Liso{t) exists 
when studying time-resolved spectra of individual bursts is a strong 
argument in favor of the reality of the spectral energy correlations, 
(Ghirlanda, Nava& Ghisellini 2010; Ghirlanda et al. 201 1) and mo- 
tivates the search for the underlying process generating them. Even 
if several ideas have been already discussed in the literature, there 
is no general consensus yet, and a step forward towards a better 
understanding both of the spectral energy correlations and the un- 
derlying radiation process of the prompt emission of GRBs is to 
discover what are the typical energetics, peak frequencies and peak 
luminosities in the comoving frame. 

The physical model of GRBs requires that the plasma emitting 
7-rays should be moving relativistically with a bulk Lorentz factor 
Fo much larger than unity. The high photon densities and the short 
timescale variability of the prompt emission imply that GRBs are 
optically thick to pair production which, in turn, would lead to a 
strong suppression of the emitted flux, contrary to what observed. 
The solution of this compactness problem requires that GRBs are 
relativistic sources. From this argument lower limits Fo ^ 100 are 
usually derived (Lithwick & Sari, 2001). The first observational 
evidences supporting this scenario were found in the radio band 
where the ceasing of the radio flux scintillation (few weeks after 
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the explosion as in GRB 970508; Frail et al. 1997), allowed to esti- 
mate r of a few. This value corresponds to the late afterglow phase, 
when the fireball is decelerated almost completely by the interstel- 
lar medium and is characterized by a much smaller bulk Lorentz 
factor than the typical To of the prompt phase. 

Large Lorentz factors imply strong beaming of the radiation 
we see. We are used to consider GRB intrinsic properties (Speak, 
Siso, iiso) for the bursts with measured redshifts, but still an im- 
portant correction should be applied. Our aim is to study the dis- 
tributions of Speak, Eiso, Lisa and the spectral-energy correlations 
(Speak — Siso and Speak — Siso) in the comoving frame, account- 
ing for the To factor. The estimate of To is possible by measuring 
the peak of the afterglow (Sari & Piran 1999) and has been suc- 
cessfully applied in some cases (e.g. Molinari et al. 2007, Gruber 
et al., 2011) and more extensively recently by Liang et al. (2010) in 
the optical and X-ray band. Other methods allow to set lower limits 
(Abdo et al. 2009; Ackerman et al. 2010; Abdo et al. 2009a) mainly 
by applying the compactness argument to the high energy emission 
recently detected in few GRBs at GeV energies by the Fermi satel- 
lite (see Zou, Fan & Piran 2011; Zhao, Li & Bai 2011; Hascoet et 
al. 201 1 for more updated calculation on these lower limits on To). 
Conversely, upper limits (Zou Sl Piran 2010) can be derived by re- 
quiring that the forward shock emission of the afterglow does not 
appear in the MeV energy band. 

The paper is organized as follows: in § 2 we discuss the rel- 
ativistic corrections that allow us to derive the comoving frame 
Speak, ^i'soand I/Jge, from the rest frame Speak, Siso, L\so', in § 3 
and § 4 we derive a general formula for the estimate of To from the 
measurement of the time of the peak of the afterglow emission; in 
§ 5 we present our sample of GRBs and in § 6 our results which are 
finally discussed in § 7. Throughout the paper we assume a standard 
cosmology with h — Qa = 0.7 and Jim = 0.3. 



spectrum, or can be the spectral peak energy of a given time in- 
terval. In this paper we will use the time-integrated Speak ~ 
^pcak(l + z). The received fluence dT/d6 (i.e. the flux inte- 
grated in time) from each annulus of same viewing angle 6 is 
dT/d9 oc sin^^'"*. For ^ — > the Doppler factor is maximum, 
but the solid angle vanishes, while for 9 > l/F the solid angle is 
large, but 5 is small. Therefore there will be a specific angle 9 for 
which dJ-/d6 is maximum. This is given by 



cos e = /? + 
At this angle the beaming factor is 

3 

We then set S^eak = Speak/(5r/3). 

he fit 
Lgqui 



5r2 



(3) 



(4) 



Siso — This is proportional to the fluence T, and the relation be- 
tween the observed and comoving quantity is 



(5) 



S;',„ r j; sin 9d9 

We then set E^^ = Siso/F. 
Siso — This is proportional to the flux S, so the ratio Liso/SLo is 



Siso_ s _j:s^ 



SJso F' Jo sin 9d9 3 



4 2 



(6) 



We then set L(so = Siso/(4r^/3) (in agreement with Wijers & 
Galama 1999). 



2 FROM THE REST TO THE COMOVING FRAME 

In this section we derive the Lorentz transformations to pass from 
rest frame quantities to the same quantities in the comoving frame. 
This is not trivial, since, differently from the analog case of blazars, 
the emitting region is not a blob with a mono-directional velocity, 
but a fireball with a radial distribution of velocities. Therefore, an 
observer located on axis receives photons from a range of viewing 
angles, complicating the transformations from rest frame to comov- 
ing quantities. We are interested to three observables: the peak en- 
ergy Speak, the isotropic equivalent energy Siso and the isotropic 
equivalent peak luminosity Siso. Dealing with isotropic equivalent 
quantities, we can assume that the emitting region is a spherical 
shell with velocities directed radially. We also assume that the co- 
moving frame bolometric intensity /' is isotropic. We then adopt 
the usual relation between observed (/) and comoving (/') bolo- 
metric intensity: 



S = 



(1) 



r{l- l3cos9) 

where 5 is the Doppler factor and 9 is the angle between the velocity 
vector and the line of sight. The received flux is 



F = 27r/' 



5'' sin ed9 



(2) 



Since the fluence J^ is a time-integrated quantity we have T oc 
J^ (5^ sin 9d9, i.e. one power of 5 less. 

Speak — This quantity can be derived from the time-integrated 



3 ESTIMATE OF THE BULK LORENTZ FACTOR Fo 

In the thin-shell regime (i.e. for Tgo < ipeak.obs, condition satis- 
fied for almost all bursts in our sample) the standard afterglow the- 
ory predicts that the peak of the bolometric afterglow light curve 
corresponds to the start of the fireball deceleration. The decelera- 
tion radius is commonly defined as the radius at which the swept up 
matter m(rdec) is smaller by a factor Fo than the initial shell's rest 
mass M() = So/(Foc^). Usually, the deceleration time idee is es- 
timated as tdcc = rdoe/(2cFo) (Sari & Piran 1999). This relation 
is approximate, since it does not consider that the Lorentz factor 
is decreasing (e.g. Bianco & Ruffini 2005). Some authors consider 
this relation to estimate Fo from the peak time of the afterglow light 
curve (Sari & Piran 1999; Sari 1997), while other authors consider 
that idee = rdee/(2cFdee), whcrc approximately Fo ~ 2r(rdce) 
(Molinari et al. 2007). 

We propose here a detailed and general calculation of Fo 
which extends the estimate to the generic case of a circumburst 
density profile described by n = nor^". We use the shape of 
the light curve in two different power-law regimes: the coasting 
phase when r <C rdoe and F(r) — Fq, and the deceleration phase 
when Tdee '^ r <^ tnr (where rNR marks the start of the non- 
relativistic regime). During the deceleration regime the evolution of 
the Lorentz factor is described by the self-similar solution found by 
Blandford & McKee (1976): 



(17 - 4s)So 



(7) 



(12-4s)m(r)c2 
The relation between the radius and the observed time is obtained 
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by integrating the differential equation dr — 2cV^{r)dt and by 
considering the exact evolution of T with r. From Eq.[6] 






(8) 



where the dissipated comoving energy -Ediss is given by (Panaitescu 
& Kumar 2000): 



iSks = (r - l)m(r)c2 



(9) 



Only a fraction Se of the dissipated energy is radiated. We assume 
that this quantity is small and does not affect the dynamics of the 
fireball (adiabatic regime). Eq.[8]holds until the emission process is 
efficient (fast cooling regime). 

During the coasting phase F = To 2> 1 and the luminosity 
(denoted by I/iso,i) is: 



T _, 4 3 2 dm(r) 



£e-FoC 47rr^ 'norrip (10) 



Since in this phase the Lorentz factor is constant and equal to Fo 
the relation between the fireball radius and the observed time is 

r = 2ciFo 

As a function of time, the luminosity is: 



r 4 (4_ji) (5-s)-n8-2s,2- 

Liso.i = ee-2 'imomp& 'Fg t 



(11) 



For a homogeneous density medium (s = 0) the light curve rises 
as t^ . The luminosity is instead constant when s — 2, which corre- 
sponds to the stellar wind density profile. 

To derive the luminosity during the deceleration phase we start 
again from Eq.[8]and Eq. |9] However, in this case F is decreasing 
according to Eq.[7](but still F 3> 1). We derive: 



r 4 2 2 

-Liso,2 — '^e-^1 C 



„dm(r) , . dr 



dt' ' ' dt' 

The first term of the sum in square brackets can be written as 

dm{r)dr_ „^ m(r) ^2„ 

dr dt' '^■^ '> r ^ " 
The second term of the sum becomes 
dF dr 
IrdP 
During the deceleration 



(12) 



m{r)- 



3-sm{r) ^2^ 



* " 2c / F2 



2(4 - s)cF2 

where we have used F(r) given in Eq.|7] 

For Fo S> 1 the initial energy content of the fireball Eq = 
Sk.iso + Mqc^ ~ i5k,iso, where Sk.iso is the isotropic kinetic en- 
ergy powering the expansion of the fireball in the ISM during the 
afterglow phase. If the radiative efficiency rj of the prompt phase is 
small, -Ek.iso can be estimated from the energetics of the prompt as 
-Ek.iso = Eiso/n- We obtain: 



Li, 



_ 4 2 2 (3-g)m(r) 
" ^^3 4(4 -s)i 

4(17-4£)(3-s)^ 1 



(13) 



3 4(12-4s)(4-s)j7 

The peak time of the light curve is the time when the coasting 
phase ends and the deceleration phase starts and can be estimated 
by setting Z/iao,i(ipoak) = -'^iso,2(tpoak): 



Cpcak — 



(17-4s)(3-s)£:i,o 



2^-^itnompC'~^^ 



''77(12-4s)(4-s)rg-2'' 
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(14) 



and inverting this relation to obtain the initial Lorentz factor as a 
function of the peak time: 



(17-4s)(3-s)Eiso 



26-''7rnompc5-»77(12 - 4s)(4 - s)t 



peak 



(15) 



where tpoak is the peak of the afterglow light curve in the source 
rest frame, i.e. fpoak = ipcak,obs/(l + z), and it will be indicated 
as tp.z hereafter. 

While a wind density profile (hereafter W: wind interstellar 
medium) is expected from a massive star progenitor that under- 
goes strong wind mass losses during the final stages of its life 
(Chevalier & Li 1999), it is not possible at the present stage to pre- 
fer the W to the homogeneous interstellar medium case (H, here- 
after). We already showed (Nava et al. 2006) that the collimation 
corrected iSpeak — E-y correlation (so called "Ghirlanda" correla- 
tion; Ghirlanda, Ghisellini & Lazzati 2004) has a smaller scatter 
and a linear slope when computed under the assumption of the W 
compared to the H case. It is, therefore, important to compare the 
estimates of Fo and of the comoving frame energetics in these two 
possible scenarios. The most extensive study of Liang et al. (2010) 
estimated Fo mostly from the peak of the afterglow light curve in 
the optical band and in few cases from a peak in the the X-ray band. 
They considered only the H case and found a strong correlation be- 
tween Fo and the ORB isotropic equivalent energy -Eiao. 

Eq. [TT] predicts that the afterglow light curve is flat in the 
coasting phase, with no peaks in the W density case (s = 2). 
However, this equation neglects pre-acceleration of the circum- 
burst matter due to the prompt emission itself, that can have im- 
portant consequences, as we discuss below. 



4 HOMOGENEOUS OR WIND DENSITY PROFILE? 

In the following we will find the initial bulk Lorentz factor Fo for 
bursts showing a peak in their early afterglow light curve. In the 
simple case of an homogeneous circumburst density, we expect that 
the afterglow luminosity I/aft on t^T^ , and therefore Laft oc t^ 
when r = Fo = constant (Eq. [TTt . It can be questioned if, in the 
case of a wind density profile, such a peak occurs, or if the initial 
light curve is flat (i.e. oc t"), as suggested by Eq.llllwhen s = 2. 

The derivation leading to Eq.[TT]assumes that the circumburst 
medium is at rest when the fireball impacts through it (i.e. it is an 
external shock). Instead, since the electrons in the vicinity of the 
burst scatter part of the prompt emission of the burst itself, some 
radial momentum has to be transferred to the medium (as suggested 
by Beloborodov 2002). If the velocity acquired by the circumburst 
matter becomes relativistic, then the fireball will produce an inter- 
nal shock when passing through the medium, with a reduced effi- 
ciency. 

To illustrate this point, let consider an electron at some dis- 
tance r from the burst, scattering photons of the prompt emission 
of energy -Bpoak = xnicC^. In the Thomson limit of the scattering 
process, this electron will scatter a number r of prompt photons 
given by: 



r — aTU-yAr 



CTiisoCiburst 

'inr^cxmcc'^ 



O'tEiso 

ATvr^xmcC^ 



(16) 



To evaluate the distance r up to which this process can be rel- 
evant, consider at what distance the electrons make a number 
r ~ {mp/me)/x scatterings, namely the distance at which the 
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GRB 


z 


-E-peak 

keV 


Siso 

erg 


erg/s 


s 


Th 


Tw 


Ref 


990123 


1.60 


2031±161 


(2.39±0.28)E54 


(3.53±1.23)E53 


18 


312 


182 


2 


030226 


1.986 


290± 63 


(6.7±1.2)E52 


(8.52±2.23)E51 


4340 


26 


19 


5 


050820A 


2.612 


1325±277 


(9.75±0.77)E53 


(91±6.8)E51 


108.17±4.62 


142 


93 


1 


050922C 


2.198 


417±118 


(4.53±0.78)E52 


(190±2.3)E51 


42 


138 


55 


2 


060210 


3.91 


575±186 


(4.15±0.57)E53 


(59.5±8.0)E51 


97 


133 


77 


2 


060418 


1.489 


572±114 


(1.28±0.10)E53 


(18.9±1.59)E51 


60.73 ±0.82 


137 


65 


1 


060605 


3.78 


490±251 


(2.83±0.45)E52 


(9.5±1.5)E51 


83.14 ±2.7 


101 


41 


1 


060607A 


3.082 


575±200 


(10.9±1.55)E52 


(20±2.7)E51 


42.89 ±0.62 


153 


68 


1 


060904B 


0.703 


135±41 


(36.4±7.43)E50 


(7.38±1.4)E50 


271.91±33.75 


50 


18 


1 


061007 


1.261 


902±43 


(8.82±0.98)E53 


(17.4±2.45E52 


34.62 ±0.18 


215 


121 


1 


061121 


1.314 


1289±153 


(2.61±0.3)E53 


(141±1.5)E51 


250 


88 


54 


4 


070110 


2.352 


370±170 


(5.5±1.5)E52 


(45.1±7.52)E50 


350 


64 


34 


4 


071010B 


0.947 


101±23 


(2.12±0.36)E52 


(64±0.53)E50 


67 


105 


40 


2 


0803 19C 


1.95 


1752±505 


(15±0.79)E52 


(9.5±0.12)E52 


117.38±3.22 


109 


57 


1 


080804 


2.2 


810±45 


(1.15±0.2)E53 


(2.69±0.32)E52 


40.5 


157 


70 


5 


080810 


3.35 


1488±348 


(3.91±0.37)E53 


(9.27±0.87)E52 


27.02 ±0.26 


214 


105 


1 


08 1203 A 


2.1 


1541±757 


(3.5±0.3)E53 


(28.1±1.94)E51 


118.09±0.46 


121 


70 


1 


090102 


1.547 


1148±143 


(2.2±0.26)E53 


(8.7±0.56)E52 


20.3 


221 


97 


5 


090618 


0.54 


155.5±11 


(2.53±0.25)E53 


(2.05±0.1)E52 


51.9 


158 


80 


5 


090812 


2.452 


2023 ±663 


(4.03±0.4)E53 


(95.6±9.66)E51 


17.38 


253 


118 


5 


091024 


1.092 


794±231 


(2.8±0.3)E53 


(1.0±0.22)E52 


1912 


59 


66 


6 


091029 


2.752 


230±66 


(7.4±0.74)E52 


(13.2±0.73)E51 


88 


111 


51 


5 


100621A 


0.542 


146±23.1 


(4.37±0.5)E52 


(3.16±0.24)E51 


3443 


26 


18 


5 


100728B 


2.106 


404±29 


(3.0±0.3)E52 


(18.6±1.20)E51 


16 


188 


63 


5 


100906A 


1.727 


158±16 


(3.34±0.3)E53 


(24.5±0.86)E51 


37 


186 


93 


5 


1 10205A 


2.22 


715±239 


(5.6±0.6)E53 


(2.50±0.34)E52 


311 


89 


62 


5 


110213A 


1.46 


241±13 


(6.4±0.6)E52 


(20.9±0.58)E51 


81 


113 


51 


5 


0809 16C 


4.35 


2759±120 


(5.6±0.5)E54 


(10.4±0.88)E53 


1.5 


880 


419 


3 


090510 


0.903 


4400±400 


(5.0±0.5)E52 


(1.78±0.12)E53 


0.44(315.3) 


773(66) 


175(34) 


3(7) 


090902B 


1.822 


2020±17 


(44±0.3)E53 


(58.9±0.97)E52 


3.2 


643 


327 


3 


090926A 


2.106 


907±7 


(20±0.52)E53 


(74±1.45)E52 


2.9 


605 


275 


3 



Table 1. The sample of GRBs with redshifts z, rest frame peak energy Speak, isotropic equivalent energy Ei^a and luminosity Liso (integrated in the 1 
keV-10 MeV energy range) and peak time of the optical afterglow light curve (given in the source rest frame tp,z). The Fq factors computed in the H and W 
case are reported. The GRBs shown separately at the bottom of the table are the three long GRBs (080916C, 090902B, 090926A) showing a peak of the GeV 
light curve (as detected by Fermi-hAI) which could be interpreted as afterglow emission (Ghisellini et al. 2010). The short GRB 090510 is shown with two 
entries: one corresponding to the peak of the GeV light curve and the second to the peak of the optical light curve. The last column gives the references for the 
peak time of the afterglow: (1) Liang et al. 2010, peak of the optical light curve; (2) Liang et al. 2010, references in their Tab. 6; (3) Ghisellini et al. 2010; e 
(4) Ghisellini et al., 2009; (5) GRBs added in this work (Melandri et al. 201 1); (6) Gruber et al. 201 1 ; (7) De Pasquale et al. 2009. 



electrons and their associated protons are accelerated to 7 ~ 2: 

-,1/2 



r(7 = 2) 



O-T-Eis, 



4:TvmpC^ 



1.9 X 10" _E, 



15pl/2 



(17) 



where -Biso,53 ~ 10' iSiso erg. This distance must be compared 
with the deceleration radius r^cc in the case of a wind density pro- 
file corresponding to a mass loss M and a velocity v„ of the wind: 



i{r) = 



M 



ATvr^mpVvj 



= 3.16 X 10'' 






(18) 



where M = W'-'M^sMq yr"^ and v^ = 10* cm s"^ (i.e. 10* 
km s^^) (e.g. Chevalier & Li 1999). The deceleration radius is 



^dcc — 



Eiso 



ATvmpC^rjrQ 



1.7 X 10i« ^'-:^*"--^ cm 

»7_iM_5rg_2 



(19) 



where 77 is the efficiency of conversion of the kinetic energy to ra- 
diation (I/iso ~ '7^k,iso)- Therefore it is possible to have a pre- 
acceleration of the circumburst matter up to a distance comparable 
to (but less than) the deceleration radius. In this case we expect to 



have a very early rising afterglow light curve (corresponding to rel- 
atively inefficient internal shocks between the fireball and the pre- 
accelerated circumburst medium), followed by a flat light curve and 
then a decay. 

We conclude that the absence of a flat early light curve does 
not exclude (a priori) a wind density profile. This gives us a mo- 
tivation to explore both cases (i.e. homogeneous and wind density 
profile) even if the bursts in our sample all show a peak in the af- 
terglow light curve (and thus a rising phase). 

Note that the same pre-acceleration can occur if the den- 
sity is homogeneous. In this case, again, we expect the very early 
afterglow to be less efficient than what predicted without pre- 
acceleration, leading to a rising phase even harder than t^ . 



5 THE SAMPLE 

Since we want to study the energetics, luminosities and peak en- 
ergies of GRBs in the comoving frame, our first requirement is to 
know the redshift z. Then we also need that the spectral peak en- 
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ergy E'^l;^ has been determined from the fit of the prompt emission 
spectrum. Most of these bursts have been localized by the Burst 
Alert Telescope (BAT; Barthelmy et al. 2005) on board the Swift 
satellite, but only for a few of them BAT could determine E'^^y_ 
(due to its limited energy range, 15-150 keV). Most of the -Bpelk 
were determined by the Konus-Wind satellite (Aptekar et al. 1995), 
or, since mid 2008, by the Gamma Burst Monitor (GBM; Meegan et 
al. 2009 with energy bandpass 8 keV-35 MeV) on board the Fermi 
satellite. Our sample of GRBs with z and constrained -Epcak {sad 
consequently with computed _Eiso and Liso) is updated up to May 
2011. It contains 132 GRBs with z, E'^^l^^ and Ei^o- We have Liso 
for all but one of these bursts. 

Within this sample, we searched the literature for bursts with 
evidence of the peak of the afterglow or an estimate of the To factor: 

(i) Liang et al. (2010 ~ LIO hereafter) measured the peaks in the 
optical light curves of GRBs and then estimated To for the H case. 
From LIO we collected 9 measurements of tp,z. LIO also collected 
other estimates of fp,z from the literature (their table 6) from which 
we get other 4 values of this observable. Therefore from LIO we 
collected 13 estimates of tp,z from the optical light curves; 

(ii) two GRBs, not included in the sample of LIO, that show a 
peak in their optical afterglow light curves are taken from Ghisellini 
et al. (2009); 

(iii) LIO searched for bursts with evidence of the afterglow 
peak up to December 2008. Our sample of bursts with redshifts, 
^peak and isotropic energies/luminosities extends to May 201 1 . We 
searched in the literature for fp,z of bursts after December 2008 and 
in 10 cases we could build the light curve with available published 
data (that will be presented in a forthcoming paper - Melandri et 
al. 201 1). Our systematic search of the literature resulted in other 2 
GRBs with a peak in the optical light curve. 

Our sample is thus composed of 27 GRBs with an estimate 
of ip.z obtained from their optical light curves. All these are long 
GRBs. 

The sample is presented in Tab. [T] where we show the relevant 
properties of these bursts used in the following sections. Col. 1 and 

2 show the GRB name and its redshift. Col. 3 the rest frame peak 
energy iJpcak, and Col. 4 and 5 the isotropic equivalent energy Ei^o 
and luminosity Liso, respectively. In Col. 6 it is reported the rest 
frame fp,z from which we compute the Fq factor in the H case (Col. 
7) and in the W case (Col. 8) assuming a typical density value no = 

3 cm~^ or no = 3 x 10"^*cm~^ (for the H and W respectively) 
and a typical radiative efficiency r) — 0.2. We note from Eg. 1151 
that the resulting Fq is rather insensitive to the choice of no and 77 
both in the H case [i.e. Fq cx (no?;)^^'*] and in the W case [i.e. 
Fo oc {novr^^% 

There are also four GRBs, detected by the Large Area Tele- 
scope on board Fermi at GeV energies, showing a peak in their GeV 
light curves (Ghisellini et al. 2010). The interpretation of the GeV 
emission as afterglow (Barniol Duran & Kumar 2009, Ghirlanda et 
al. 2010, Ghisellini et al. 2010) is however debated (Ackermann et 
al. 2010; Piran & Nakar 2010). Among these bursts there is also the 
short/hard GRB 090510 whose Fo is derived from the modeling of 
the GeV light curve (Ghirlanda et al. 2010a). However, this burst 
also shows a clear peak in the optical at ~300 s after the GRB onset 
(De Pasquale et al. 2009) which questions the afterglow interpreta- 
tion of the GeV emission. 

The three LAT bursts with ip,z measured from the GeV light 
curve and the short GRB 090510 are shown separately in Tab. [T] 
These events have the smallest fp,z in our sample and, therefore, the 
largest Fo values (see Tab.lTJ. This is expected since, as discussed 
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Figure 1. To distiibutions of the 31 GRBs in the case of an homogeneous 
interstellar medium (H - solid filled blue histogram) and in the case of a 
wind density profile (W - hatched histogram). The Gaussian functions show 
the fits (solid and dashed line for the H and W case, respectively) to the 
histograms of the sample of 27 GRBs with fp_z derived from the optical 
light curve. The three long and one short GRBs with tp,z measured from 
the GeV light curve are shown by the grey solid and hatched histograms, 
for the H and W case respectively, but are not included in the fits. 



in Ghisellini et al. (2010), the detection in the GeV energy range by 
LAT seems to be a characteristic of GRBs with the largest values 
of Spclk- Besides, the possible measure of f p,z in the optical range 
is limited by the time delay of the follow up of GRBs in this band, 
although several GRBs have been repointed in the optical band by 
UVOT on board Swift. In the end, there could be a selection bias 
on the bursts with a peak in the GeV energy range, coupled with 
the debated interpretation of the GeV emission as afterglow. For 
these resons, in the next sections we will present the results of the 
study of the correlations between the GRB energetics and Fo both 
including and excluding these bursts. In all our quantitative analysis 
we always excluded the short GRB 0905 10 which is only shown for 
comparison with the properties of the 27 long GRBs. 

In our sample we do not include upper limits on tp,z which 
are those bursts observed early in the optical whose light curve is 
decaying up to several days without any sign of a peak. Several of 
these cases can be found in the literature and they would provide 
lower limits on the value of Fo. However, it is hard to define an 
appropriate sample of upper limits on ip,z derived from the optical 
band because of the lack of a unique follow-up program dedicated 
to the systematic observations of GRB afterglows. 



6 RESULTS 

In this section we first show the distributions of the Fo factors com- 
puted in the H and W and show the correlation of Fo with the 
isotropic energy Eiso and luminosity Liso. Then we show how the 
distributions of Ep^ak, L^iso and Liso change when they are cor- 
rected for the Fo factor, i.e. how they appear in the comoving frame 
(-^'poak' ^-i'so' ^iso)- In dolug thls We always consider the two es- 
timates of Fo in the H and W to compare the different distribu- 
tions of the spectral parameters. Finally, we present the rest frame 
L^peak — L^iso and L'pcak " Liso Correlations (updated here with 
132 and 131 GRBs up to May 2011) and, for those bursts in our 
sample with measured Fq, we show where they cluster in these 
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Paiameter 
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132 


2.68 


0.43 
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2.81 


0.50 
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2.85 


0.35 
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132 
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0.77 
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log -Liso 


131 
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Density 


H log To 
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27 
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27 
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0.55 
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27 
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30 
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0.39 


W log To 


27 


1.82 


0.20 




30 


1.82 


0.21 


l°g-E^;oak 


27 


0.79 


0.24 




30 


0.76 


0.27 


logi^Lo 


27 


51.47 


0.43 




30 


51.54 


0.45 


log^L 


27 


48.69 


0.26 




30 


48.71 


0.23 



6.2 E\so-Tq, Liso-ro, S| 



peak" 



-Fo correlations 



Table 2. Central values and dispersions of the Gaussians fitted to the distri- 



butions of Tq, Epeak 



^d-E^Uk'-Ei. 



, and E'- , Liao and L'- For each 



quantity we report the Gaussian fits to the sample of 27 GRBs with tp_z 
measured from the optical fight curve and the sample of 30 GRBs which 
includes the three events with tp,z measured from the GeV light curve, if 
interpreted as afterglow. The short GRB 090510 has been excluded from 
this analysis. 



planes when the beaming corrections (Bpeak = i5pcak/(5r/3), 
EL = Siso/r, LL = Liso/(4rV3)) are applied. 

For all the reasons outlined in §5, in the following we consider: 

— the optical sample of 27 GRBs with measured z, E'^l^, E\so 
and Liao, whose tp,z is measured from the optical light curve. 

— the extended sample of 30 GRBs which includes the three long 
GRBs with a peak in the GeV which, if interpreted as afterglow 
emission, allows to estimate the largest Fo in our sample. 



In this section we explore the presence of correlations between the 
rest frame GRB properties (i.e. the peak energy Speak, the isotropic 
equivalent energy Eiso and luminosity Liao) and the Fo factor 

In the upper panels of Fig. |2] we show the isotropic energy 
Eiso and luminosity Lisa (open red circles and filled green squares, 
respectively) as a function of Fo in both the H and W case (left and 
right panel, respectively). In the bottom panels of Fig. [2] we show 
the peak energy Speak as a function of Fo in the H (left panel) and 
W (right panel) case. 

Fhe Spearman rank correlation coefficients and associated 
chance probabilities are reported in Tab. [3] We model the corre- 
lations with a power law: logY = m logFo + q (with Y=Eiso, 
Y=Liao or Y=Epaa,k) and list the best fit parameters in Tab. [3] We 
fit this model to the data points (shown in Fig. [2J with the bisec- 
tor method. The choice of this fitting method, instead of the least 
square Y vs. X method that minimizes the vertical distances of the 
data from the fitting line, is motivated by the large dispersion of 
the data and the absence of any physical motivation for assuming 
that Fo or instead Eiso, iiso or Speak are the independent variable 
(Isobe et al. 1990). 

In a recent work, Lv et al. (2011) derive a correlation Fo oc 
Efso'^, similar to that found in LIO. Such a flat correlation is ob- 
tained because Fo is fitted versus Eiso (or Liao)- As described 
above, the large scatter of the correlations and the lack of any phys- 
ical reason for assuming either Fo or Siao (Siso) as the indepen- 
dent variable, requires instead that these correlations are fitted with 
the bisector method. This gives different correlation slopes with re- 
spect to those reported in LIO and Lv et al. (201 1). Moreover, in our 
sample we only consider bursts with firm estimates of Speak and 
do not include those GRBs which are fitted by a simple power law 
in the BAT energy range but whose peak energy is derived through 
a Bayesian method, based on the properties of bright BATSE bursts 
(Butler et al. 2008). 

We find that there are strong correlations between the spectral 
peak energy and isotropic energy/luminosity with Fo. The slopes 
of these correlations are rather insensitive to the circumburst pro- 
file adopted in deriving Fo (H or W) and are similar for Siao and 
Siso (SisoOcFo^ and LiaoOcFo^). A roughly linear correlation exists 
between Speak and Fo: SpeakOcFo (bottom panels in Fig.[2]l. 

The dispersion of the data points around the best fit corre- 
lations (shown by the solid and dashed lines in Fig. |2j is mod- 
eled with a Gaussian and its CTae is given in Tab. [3] The less 
dispersed correlation is between the luminosity Liao and Fo(with 
o-ac = 0.07). 

We finally verified that there is no correlation between the 
GRB duration T90 and Fo (chance probability P = 0.3 and 
P — 0.7 for the H and W case) and between the redshift z and 
Fo. 



6.1 Fq distributions 

Fig. [I] shows the distributions of the Fo factors of the 27 GRBs 
of our sample (with ip,z measured from the optical light curve - 
Tab.l) computed in the H (solid histogram) and W case (hatched 
histogram), respectively. The two distributions are fitted with Gaus- 
sian functions and the central value and dispersion are reported in 
Tab.|2] The average Fq factor is ~138 in the H case and ~66 in the 
W case. In both the H and W case the distribution of Fo is broad, 
spanning nearly one decade. 



6.3 Comoving frame Spoak» ^iso» ^Lo distributions 

In Fig.[3]|4]and[5]we show the distributions of the comoving frame 
peak energy, isotropic equivalent energy and luminosity. In Fig. |3] 
we show the distributions of the peak energy: the sample of 132 
GRBs with measured redshifts and known Speak is shown with the 
dashed line and the subsample of 30 GRBs of this work for which 
we could estimate Fo is shown with the red hatched histograms. 
These distributions represent Speak, i.e. the peak energy in the rest 
frame of the sources. 

The distributions of the comoving peak energy [derived as 
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Figure 2. Top panels: Isotropic equivalent energy E^^^ (open circles) and luminosity L^^^ (filled squares) as a function of Fq, computed for the 30 GRBs in 
our sample in the H case (left panel) and W (right panel). The solid (dashed) line in both panels show the least square fit with a power law to the E^^a-^o 
(^iso-ro) correlation to the sample of 27 GRBs with peak in the optical Ught curve (open red circles and filled green squares). The three GRBs with peak in 
the GeV light curve are shown with the grey symbols but are not included in the fits shown here. The short GRB 090510 with both a peak in the GeV and a 
delayed peak in the optical (see Tab. [T) is shown by star symbols connected by the dashed (gray) line. The larger value of Fq is that derived from the peak in 
the GeV light curve. Bottom panels: Peak energy Bpealc for the H case (left panel) and W case (right panel) as a function of Fq. The solid line is the best fit 
correlation. The correlation coefficient and the slope and normalization of the best fit correlations are reported in Tab.|3] 



Correlation 


#GRBs 


P 


-^chance 


771 


9 


O-sc 




27 


0.48 


10-2 


1.92±0.40 


49.20±0.88 


0.28 




30 


0.74 


2.5 X 10-'' 


1.96±0.26 


49.11 ±0.62 


0.23 


^iso ~ Tq 


27 


0.64 


3 X 10-4 


2.15±0.34 


48.01±0.74 


0.18 




30 


0.74 


3 X 10-6 


2.04±0.22 


48.21±0.51 


0.20 


i^poak - rH 


27 


0.45 


10-2 


1.31±0.2 


0.03±0.36 


0.21 




30 


0.56 


10-3 


1.13±0.13 


0.36±0.31 


0.23 


p. pw 


27 


0.75 


4 X lO-'' 


2.36±0.36 


48.97±0.60 


0.18 




30 


0.82 


2.2 X 10-** 


2.15±0.20 


49.32±0.42 


0.10 


-^iso ~ Tq 


27 


0.76 


5 X 10-6 


2.40±0.24 


48.14±0.43 


0.07 




30 


0.82 


2.6 X 10-** 


2.19±0.16 


48.52±0.31 


0.10 


-C-pcak J- 


27 
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profile (W). The Spearman correlation coefficient p and the chance probability Pctianco ^^ reported together with the slope m and normalization q of the fit 
of the data points with a linear model The fit is done with the bisector method considering the sample of 27 GRBs with optical peak and the 30 GRBs (i.e. 
including the three long bursts with peak in the GeV). 



^pcak=^pcak/(5ro/3)] are shown by the (cyan) filled and hatched 
(purple) histograms in Fig. [3] for the H and W case, respectively, 
considering the 27 GRBs which show a peak in the optical light 
curve. Fig. [3] shows also the fits with Gaussian functions: their pa- 
rameters are reported in Tab.|2] 



There is a reduction of the dispersion of the distribution of 
the peak energy from the rest frame to the comoving one. In the 
comoving frame iJpoak clusters around ~6 keV and ~3 keV in the 



H and W case, respectively, with dispersions of nearly one decade, 
i.e. narrower than the dispersion of Speak- 

Fig. |4] shows the distribution of the isotropic energy Eiso for 
all the 132 GRBs with known z and measured Bpeak (dashed line) 
and for the 30 GRBs with an estimate of To (hatched red his- 
togram). The i5i'sQ=i?iso/ro distributions are shown with the solid 
filled (cyan) histogram and the hatched (purple) histogram for the 
H and W case. These distributions are obtained with the 27 GRBs 
with a peak in the optical light curve. The three GRBs with a peak 
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Figure 3. Peak energy distributions in the rest frame -Bpeak (dashed his- 
togram) for the sample of 132 GRBs with known redshift and constrained 
Bpeak- The hatched histogram shows the 30 GRBs of our sample for which 
we have an estimate of the peak of the afterglow and hence of To. The 
beaming corrected distribution of E' j,jjjj=i?pcak/(5ro/3) is shown by the 
solid filled (cyan) histogram in the H case and with the hatched (purple) his- 
togram in the W case. For all the distributions we also show the Gaussian 
fits whose parameters are reported in Tab.|2]- The four GRBs with a peak in 
the GeV light curve are shown with gray filled and hatched histograms. 



Figure 4. Isotropic energy distributions in the rest frame (dashed his- 
togram) for the sample of 132 GRBs with known redshift and constrained 
^°cak' '^^^ hatched histogram shows the 30 GRBs of our sample for which 
we have an estimate of the peak of the afterglow. The beaming corrected 
distribution of E[^^=E\^alT oi's, shown by the solid filled histogram and 
hatched purple histogram for the H and W case for the 27 GRBs with a 
peak in the optical light curve. The four GRBs with a peak in the GeV fight 
curve are shown for comparison with the hatched and filled gray histograms. 



in the GeV light curve are only shown for comparison (hatched and 
filled gray histogram). The distributions of E[^^ are wide. On aver- 
age the comoving frame Bi'sq^I-Sx 10^^ erg in both the H and W 
case, but there is a reduction of the dispersion of the distribution of 
Eiso from the rest {a^c ~ 0.64) to the comoving frame (ctsc = 0.43 
and (Tsc = 0.49) for the W and the H case, respectively (see Tab. 

Finally, in Fig. \5\ we show the distribution of Liso for the 
131 GRBs in the sample (dashed line), the distribution of Liao for 
the 30 GRBs with estimated To (red hatched histogram) and the 
comoving frame I/(so=Liso/(4ro/3) distribution (solid filled cyan 
and hatched purple histograms for the H and W case, respectively, 
obtained with the 27 GRBs with a peak in the optical light curve). 
Interestingly, we find a strong clustering of the comoving frame dis- 
tribution of I/iso. For the H case we find (see Tab. [2] for the values 
of the Gaussian fits) an average L[^a^ 10^* erg s~^ with a small 
dispersion (0.47 dex), while when using the To computed in the 
wind density profile (W) case we find an almost universal value of 
Ljgo~ 5 X 10*** erg s~^ with a dispersion of less than one order of 
magnitude around this value (hatched purple histogram and dashed 
purple line in Fig.[5](. 



We have 132 GRBs with measured z and known Speak and Eiso 
and 131 GRBs with measured z and Speak and Liso- We show 
the corresponding Speak — Eiso and Speak — Siso correlations 
in Fig. |6] (left and right panel respectively). The best fit correla- 
tion parameters (obtained with the bisector method) are reported in 
Tab.H We find that Speak oc S^^'^ (dashed line in Fig. [til with a 
scatter a = 0.24 (computed perpendicular to the best fitting line 
and modeled with a Gaussian function). The other correlation is 
Speak oc Li^o with a slightly larger scatter a = 0.3. The 1, 2 and 
3(T dispersion of the correlations are shown with the shaded stripes. 
Fig. [6] also shows the comoving frame E'^^^^ and E{^^ (left 
panel) and Speak and L[g^ (right panel) for the 30 GRBs of our 
sample with an estimate of To in the H case. The 27 GRBs with a 
peak in the optical are shown with the cyan filled squares in Fig|6] 
while the three long GRBs with a peak in the GeV light curve are 
shown with the filled gray squares. Fig. |7] show the same correla- 
tions (Speak — Siso and Speak ~ ^iso in the left and right panels 
respectively) for the W case. We note that in both the H and W 
cases there is a clustering of the points around typical values of 
Speak> -Eiso and L[^a- Tab.|4]reports the correlation analysis among 
the comoving frame quantities. 



6.4 Comoving frame Speak 
correlations 



EU and E' 



peak 



T' 



Here we show the effect of correcting the spectral energy correla- 
tions Speak — Siso and Speak — iiso for the bulk Lorentz factors 
To . These correlations were originally found with a dozen of GRBs 
(Amati et al. 2002 and Yonetoku et al. 2004 for the Speak - Siso 
and Speak — Siso correlations respectively) and since then updated 
with newly discovered GRBs with measured redshifts z and well 
constrained spectral peak energies Speak- In this work we have up- 
dated the sample of GRBs with all these observables to May 201 1. 



7 DISCUSSION AND CONCLUSIONS 

We have considered all bursts with measured Speak and known red- 
shift up to May 2011 (132 GRBs). Among these we have searched 
in the literature for any indication of the peak of the afterglow light 
curve tp,z suitable to estimate the initial bulk Lorentz factor Fq. 
Our sample of bursts is composed by 27 GRBs with a clear evi- 
dence of ip,z in the optical light curve. We have derived the peak 
energy Speak, 'h^ isotropic energy E^^^ and the isotropic peak lu- 
minosity Ljso in the comoving frame. To this aim we have derived 
the general formula for the computation of To (§.3) considering 
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Figure 5. Isotropic luminosity distributions in the rest frame (dashed histogram) for the sample of 131 GRBs with known redshift and constrained E°^^^^. 
The hatched histogram shows the 30 GRBs of our sample for which we have an estimate of the peak of the afterglow. The beaming corrected distribution of 
L' is shown by the solid filled histogram and hatched purple histogram for the H and W case for the 27 GRBs with a peak in the optical light curve. The 
four bursts with a peak in the GeV light curve are shown for comparison with the hatched and filled grey histograms. 
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Table 4. Results of the fit of the Epeak — ^iso and Epeak — ^iso correlations updated in this paper to May 2011. The Spearman correlation coefficient p 
and the chance probability Pchanco is given with the slope m and normalization q of the least square fits. 



© 0000 RAS, MNRAS 000, 000-000 



10 G. Ghirlanda et al. 



10 




10- 



10 



10 10- 

E. [erg] 



10- 



10^ 



10^' 



10^' 



10^' 
L [erg/s] 



10 



W 



Figure 6. Homogeneous interstellar medium - H. Left: -Epeak ~ ^iso correlation in the rest frame (crosses and red circles) for 132 GRBs with z and fitted 
E-peak updated to May 201 1. Right: -Bpoak ^ iiso correlation with 131 GRBs. In both panels the best fit correlation is shown by the dashed line and its 1, 2, 
3(7 scatter is shown by the shaded region. The comoving frame E' ^^y^ and E'^^^ (left) and E' ^^^ and LJ^^ (right) of 30 GRBs (red open circles [left panel] 
and green open circles [light panel]) in our sample (Tab.[T) with an estimate of the Fq factor are shown with the filled cyan square symbols (27 events with 
tp^z in the optical light curve) or grey filled square (the three long GRBs with a peak in the GeV light curve). The short GRB 090510 is also shown with a star 
symbol and the low luminosity GRB 060218 (with ro~5 [Ghisellini et al. 2006]) is shown with an open circle. 
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Figure 7. Wind interstellar medium - W. Same as Fig.|6] 



two possible scenarios: a uniform interstellar medium density pro- 
file (n =const, H) or a wind density profile (n oc r~^, W). 

For the wind case the To -distribution (Fig. [T] and Tab. |2j is 
shifted at somewhat smaller values ((Fq) ~ 66) than the same dis- 
tribution for the homogeneous density case ((Fo) ~ 138). The dis- 
tribution of iJpcak is relatively narrow and centered around ~6 keV 
or ~ 3 keV for the W and H case (Fig. [5] and Tab. [2j- The distri- 
bution of LJso (Fig. |5]l clusters, especially for the wind case, in a 
very narrow range (much less than a decade), around 5 x 10** erg 



s^^, while the distribution of El^^ (Fig.|4j is broader and centered 
at 3 X 10^^ erg. Eiso and Liso correlate with Fo, (ocFo^'^ both for 
the wind and the homogeneous case) and the correlation is stronger 
(with a scatter a = 0.07) for the wind case. Finally, the duration of 
the burst, as expected, does not correlate with Fq. 

The correlations that we have found are strong despite they 
are defined with a still small number of GRBs. We expect that with 
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determined spectral properties (i.e. -Bpeak, ^iso and Liso) the slope 
and normalization of these correlations might change. 

For comparison we also considered four GRBs with a peak in 
the GeV light curve. If the GeV emission is interpreted as afterglow 
(Barniol-Duran & Kumar 2009; Ghirlanda et al. 2010; Ghisellini 
et al. 2010) the measure of ip,z at early times in the GeV range 
allows us to estimate their To, that are consistent with the correla- 
tions found using only the bursts with fp,z observed in the optical. 
Although not a proof, this is a hint in favour of the afterglow origin 
of the GeV emission. 

These results are schematically summarized in the first col- 
umn of Tab. |5] The second column of the same table reports some 
immediate implications of these results. Since -Epeak '^ ^pcakTo is 
contained in a narrow range, all bursts emit their radiation at a char- 
acteristic frequency in their comoving frame, irrespective of their 
bulk Lorentz factor. Furthermore, we can assume that -Bpeak ocFq, 
and this, together with the quadratic dependence on Fo of iJiso and 
Liso, yields the "Amati" and the "Yonetoku" relations. They are 
the result of a different To-factors. Indeed, at the extremes of the 
Speak — EiBo and Bpeak — Liso Correlations we find GRB 060218 
which has the lowest ro~ 5 (inferred from its X-ray and optical 
properties - Ghisellini, Ghirlanda & Tavecchio 2007), while at the 
upper end (corresponding to the largest peak energies and isotropic 
energetics and luminosities) there is GRB 0809 16C which has the 
largest ro=880. The fact that the -Epcak— L^iso and -Epcak— Liso cor- 
relations could be a sequence of Fo factors has been also proposed 
by Dado, Dar & De Rujula (2007) based on different assumptions. 

If all bursts had the same jet opening angle, then L'^ = OfL[^^, 
and the (logarithmic) width of the L(go distribution would be the 
same of the (more fundamental) Li.^ distribution. On the other hand, 
we have some hints that very energetic and luminous GRBs tend to 
have narrower opening angles (e.g. Firmani et al. 2005). It is this 
property that makes the collimation corrected E-, and L^ quantities 
to correlate with JSpcak in a different way (i.e. different slope) than 
in the Amati and Yonetoku relation (Ghirlanda et al. 2004; Nava et 
al. 2006). 

We are then led to propose the following ansatz: the opening 
angle of the jet inversely correlates with the bulk Lorentz factor 
6j oc Fo^". There are too few GRBs in our sample with measured 
6j to find a reasonable value for the exponent a, but it is never- 
theless instructive to explore the case a — 1/2, leading to 9fro = 
constant. If we assume this relation we find, for the collimation 
corrected E^ : 
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Figure 8. Jet opening angle as a function of Fq for a H (stars) and for 
a W (squares). Empty symbols show the jet angles estimated by assuming 
the consistency of our sample with the -Bpcak^^7 relation. Filled symbols 
refer to the bursts of our sample for which the jet opening angle has been 
calculated from the measured jet break time of the optical light curves. The 
two lines (dashed for the H case and dot-dashed for the W case) show the 
powerlaw fit of the data points considering 9jot vs Fq and Fq vs 9j^i . The 
gray symbols show the three long bursts with a peak in the GeV light curve 
that, if interpreted as afterglow emission, allows us to estimate Fq . 



amount of energy at the same peak frequency, irrespective of the 
bulk Lorentz factor. For larger Fo the emitting time in the comoving 
frame is longer (by a factor Fo if the observed Tgo is the same), so 
the comoving luminosity is smaller. But since the jet opening angle 
is also smaller (for larger Fo), the isotropic equivalent luminosity 
turns out to be the same. These consequences are listed in the third 
column of Tab. [5] 

Interestingly, we note that the general formula for the estimate 
of the jet opening angle 



^j ,obs 
l + Z 






(24) 



with s = for the homogeneous case and s = 2 for the wind case, 
can be combined with Eq.[T5]to give: 



Ery = 6j Eiso OC Fo OC Bpeak 



(20) 



This is the "Ghirlanda" relation in the wind case (Nava et al. 2006). 
Similarly, for the collimation corrected luminosity (Ghirlanda, 
Ghisellini & Firmani 2006): 



Lry = 6j Liso OC Fo (X Speak 



(21) 



Another important consequence of our ansatz is that, in the comov- 
ing frame, the collimation corrected energetic E'^ becomes con- 
stant: 



E^, — 6\ ^° — constant 

1 



(22) 



This allows to "re-intepret" the constancy of LJso as a consequence 
of the constant E': 



L', 



so fji/ n2 



E^ 



constant 



(23) 



T^o^j rooefFo 
In other words, in the comoving frame, the burst emits the same 



jFo OC 



j.obs 



(25) 



^ ^p,obs , 

The product 6]To then depends only on two observables, i.e. the 
time of the peak of the afterglow f p,obs and the time of the jet break 
f j.obs, and it is independent from the redshift z and the energetic 
Siso as well as from the density profile normalization tiq and radia- 
tive efficiency r;. If also the product 0-^ Fo =const, then we can de- 

3-3 3-a 

rive both 6*] oc (tp,obs/ij,obs) **-== and Fo oc (tj,obsAp,obs)''-= ■ If 
the ansatz Sf Fo = const will prove to be true, then by simply mea- 
suring the peak time and the jet break time of the afterglow light 
curve we could estimate both 6\ and Fo for any GRB. 

In our sample, only for 4 bursts we can estimate the jet open- 
ing angle from the measure of the jet break time of the optical light 
curve. Their small number does not make possible to directly test 
the existence of a relation between Fo and 6'j. However, an esti- 
mate of the jet opening angle can be possible by assuming that 
all bursts in our sample are consistent with the "Ghirlanda" re- 
lation. Fig. [8] shows the estimated 6'j as a function of Fo. Stars 
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Our results 


Implications 




If 6*2 r ~const 


S^cak ~ '^onst 


-Epoak <x r 






Ei.o <x r2 


Si- « S2^^, 
Liso 0^ -Epcak 




E^, = efEi,, oc r oc Epeak 


-Liso oc i 




L^ = 6»?Liso oc r oc Speak 


Tgo not /(r) 


T^o<xr 




E!^ ~ const 


■^Lo ~ const 


KolL[,o « 7^90 


ocr 


l; ~ Ei^/T^, ~ 1/r 



Table 5. Schematic summary of our results and their implications for the case of a wind density profile. We have assumed that both -Eiso and Lis 
r^, instead of r2-2. 
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Figure 9. Distribution of 6j To in the H and W case (blue and purple his- 
tograms) estimated by assuming the -Bpcak^^7 relation in the H (Ghirlanda 
et al. 2004) or W (Nava et al. 2006) case. The hatched histograms show the 
few GRBs in our samples for which 9j has been calculated from the mea- 
sured jet break time in the optical light curve. 



(squares) refers to angles derived under the assumption of a H 
(W). To estimate the jet opening angles we considered the most 
updated "Ghirlanda" correlation, which comprises 29 GRBs with 
measured jet break time (Ghirlanda et al. 2006). For the homo- 
geneous density profile the relation has the form log iSpcak = 
-32.81 + 0.70 log £^, while in the case of a W the relation be- 
comes log Speak = —50.08 + 1.04 log _E-y. Given the large scatter 
of the data points in Fig.jS] we fitted both Q^ versus To and To ver- 
sus Sj: we obtain 6j oc F^"'^'^ and To ex: S"^"^^ for the H case 
(dashed lines in Y\%^ and 6lj oc F(7"-^^ and Ft, oc Qr^'^'^ for the 
W case (dot-dashed line in Fig. [§}. We conclude that our ansatz 

1/2 

6j OC Fq is consistent with, but not proven by, this analysis. 

An interesting exercise is to estimate the product ^jFq. From 
the observational point of view Q-^V^X at the end of the prompt 
phase, so that the decrease of F in the afterglow phase, due to the 
interaction of the GRB fireball with the interstellar medium, gives 
rise to a jet break when 6j F~ 1 . 

Some numerical simulations (Komissarov et al., 2009) of jet 
acceleration have shown that a magnetic dominated jet confined by 
an external medium should have SjFq^ 1. This value is inconsis- 
tent with typical values of Sj and Fo: in the case of an homogeneous 
wind density profile the typical Sj ~ 0.1 radiants (Ghirlanda et al. 
2007) while in the case of a wind density profile Sj ~ 0.07 radiants. 
Combining these values with the average values of Fo estimated in 
this paper (Tab.[B we find 6ljro~ 14 (5) for the H (W) case. 



These are approximate values: the sample of GRBs with mea- 
sured Sj (Ghirlanda et al. 2007) contains only 4 bursts of the sample 
of events of the present paper with estimated Fq. However, though 
somehow speculative, we can derive Sj for the 32 GRBs of our sam- 
ple assuming the Speak — S7 correlation in the H case (Ghirlanda 
et al. 2004) or in the W (Nava et al. 2006). In Fig.|9]we show the 
distributions of the product 6'jFo in the H case (blue histogram) 
and in the W case (purple histogram). We note that both are cen- 
tered around typical values of 20 and 6 (for the H and W case, 
respectively). These values are in good agreement with the results 
of recent simulations of (i) a magnetized jet confined by the stellar 
material that freely expands when it breaks out the star (Komis- 
sarov, Vlahakis & Koenigl 2010) or (ii) a magnetized unconfined 
split-monopole jet (Tchekhovskoy, McKinney & Narayan 2009; 
Tchekhovskoy, Narayan & McKinney 2010). A possible test of 
these two scenarios could be short GRBs where the absence of the 
progenitor star would prefer model (ii) for the jet acceleration. In 
our sample only the short/hard GRB 090510 is present. No jet break 
was observed for this event and in general we do not yet know if 
short GRBs follow the same Speak — S-, correlation of long ones. 
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